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Abstract 

Texture zeroes in neutrino mass matrix Mjy may give us hints about the symmetries involved in 
neutrino mass generation. We examine the viability of such texture zeroes in a model independent 
way through a bottom-up approach. Using constraints from the neutrino oscillation data, we 
develop an analytic framework that can identify these symmetries and quantify deviations from 
them. We analyze the textures of Mjj as well as those of Mm, the mass matrix of heavy Majorana 
neutrinos in the context of Type-I seesaw. We point out how the viability of textures depends 
on the absolute neutrino mass scale, the neutrino mass ordering and the mixing angle ^13. We 
also examine the compatibility of discrete flavor symmetries like /i-r exchange and 5*3 permutation 
with the current data. We show that the — r exchange symmetry for My can be satisfied for any 
value of the absolute neutrino mass, but for M^ to satisfy the 5*3 symmetry, neutrino masses have 
to be quasi-degenerate. On the other hand, both these symmetries are currently allowed for Mm 
for all values of absolute neutrino mass and both mass orderings. 
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I. INTRODUCTION 



The current low energy neutrino oscillation data [l| indicate that all three of the physical 
left-handed neutrinos have different masses and they mix among themselves. If the neutrinos 
are Majorana, the neutrino mass matrix My in the flavor basis is diagonalized by the unitary 
Pontecorvo-Maki-Nakagawa-Sakata matrix f/pMNS through 



/Kfdiag _ 7-rt T\J Tj* 



I.e. 



[/pMNS '"^^ ^PMNS 



(1) 



For Majorana neutrinos, ?7pmns is given by 
/ 



PMNS 



C12C13 



~'5l2C23 ~ Ci2S23Si3e ' 



S12C13 sise 



iS 



\ 



■S12S23 ~ Ci2C23Si3e 



C12C23 ~ -5125233136 S23C13 



"C12S23 ~ ■Sl2C23'Sl3e C23C13 



(2) 



where Cij and Sij stand for cos 9ij and sin 9ij respectively. Here f/^ = diag(e*'^S e"^^, 1), with 
the Majorana phases 0i and 02 defined in such a way that the diagonal elements of M^^^^ 
are given by 

M,^^^s = diag(mi,m2,m3). (3) 

Here mj {i = 1, 2, 3) correspond to the neutrino masses, which are chosen to be real and 
positive. The Dirac phase 6 accounts for the charge parity (CP) violation in the lepton 
number conserving processes. The phase matrix = diag(e*^'=, e*^'', e*^^) consists of the 
three "flavor phases"^ Xa that correspond to the multiplication of a neutrino flavor eigenstate 
Ua by e*^". Note that once the mixing angles 6ij have been defined to be in the first quadrant, 
the Dirac phase 6 can take values in [0, 27r) and the phases Xa, 4'i can take values between 
[0, vr). All the angles and phases are then uniquely defined. 

A global analysis of the current neutrino oscillation data at 3a CL. yields [l| 



0.25 < sin^^is < 0.37, 



0.36 < sin2^23 < 0.67 



sin2 0i3 < 0.056 . 



(4) 



^ These phases are often referred to in the Uterature as "unphysical phases" . Though these phases have 
no relevance for the low energy neutrino phenomenology and cannot be determined through low energy 
measurements, their values may be predictable within the context of specific models with new physics at 
the high scale. 
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While the absolute mass scale of the neutrinos is not yet fixed, the two mass- squared differ- 
ences have already been determined to a good degree of accuracy: 



Am| = ml-ml = (7.06 ■ • ■ 8.34) x 10"^ eV^ , 

2 



Aml^ ^ml- ^!!^1±J^^ = ±(2.07 ■ ■ ■ 2.75) x 10"^ eV^ . (5) 

It is not known whether the neutrino mass ordering is normal (mi < m2 < m^) or inverted 
(m3 < nil < ^2)- The Dirac phase 5 and Majorana phases 0i^2 are completely unknown. 
The absolute values of neutrino masses cannot be probed by oscillation experiments, the 
direct limit on the neutrino mass scale mo is obtained by the tritium beta decay experiments 
^ as mo < 2.2 eV. The most stringent constraint on mo however comes from cosmology: 
the WMAP data implies ^ 

J]m, <leV. (6) 

In this paper, we shall take the upper bound on each neutrino mass conservatively to be 
rrii < 0.5 eV. 

Given the absolute values of the neutrino masses and the complete matrix ?7pmns, the 
neutrino mass matrix Mi, in the fiavor basis can be reconstructed through Eq. ([T]). The 
structure of this matrix may reveal the presence of fiavor symmetries in the neutrino sector. 
In this paper, we consider multiple texture zeroes of M^, jg, Q, Q] as well as symmetries like 



the II — T exchange |9| and 5*3 permutation lOj, which predict certain relations between the 
elements of My. 

The symmetry-based relations among the elements of M^ and texture zeroes of My have 

n 

been explored earlier mainly by adopting a top-down approach In this approach an 

appropriate symmetry is imposed on the neutrino mass matrix, which in turn gives a pre- 
diction for the neutrino mixing parameters that can then be checked against the available 
data. We take the bottom-up approach, starting with our current knowledge about neutrino 
masses and mixings, and checking if a certain texture zero combination or a symmetry-based 
relation is allowed. This allows us to test in a model independent manner the symmetries 
present in neutrino mass generation mechanisms. It also enables us to determine which 
future measurements can act as tests of these symmetries. 

In the present approach the elements of My are expressed as functions of the absolute 
neutrino masses, the mixing angles as well as the Dirac, Majorana and fiavor phases. Our 
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current complete ignorance about these phases allows a lot of freedom for the elements of Mi, 
in spite of the relatively well measured values of the masses and mixing angles. Even with 
this freedom, some of the texture zero combinations and symmetries are clearly forbidden, as 



has been numerically verified 12j. We develop an analytical treatment, using perturbative 
expansion in appropriate small parameters, and demonstrate the analytical rationale behind 
the ruling out of some of these relations. This also leads us to the result that the additional 
knowledge of the absolute mass scale of the neutrinos and the mixing angle 613 will be crucial 
in testing for these relations in near future. 

The seesaw mechanism |l3| is one of the most favored and explored mechanisms for 
neutrino mass generation, which gives rise to light Majorana neutrinos that can satisfy the 
low energy neutrino oscillation data, as well as to heavy Majorana neutrinos that may play 



an important role in leptogenesis IJ] . If the neutrino masses are generated from a Type-I 
seesaw mechanism where three singlet heavy Majorana neutrinos are added to the Standard 
Model (SM), then we have the effective neutrino mass matrix 

= -mDM-\mDf , (7) 

where mz) is the Dirac mass matrix of neutrinos, and Mm is the Majorana mass matrix for 
the right-handed heavy Majorana neutrinos. If the heavy Majorana neutrinos are written in 
a basis where m^) is diagonal, the texture zeroes as well as symmetry relations between ele- 
ments of Mm can be related to those of the inverse neutrino matrix M~^ in a straightforward 
manner 15|] . The same analytical treatment developed for can then be extended to test 
the symmetry relations for Mm in this basis. We perform this analysis, with a particular 
emphasis on the dependence of these relations on the absolute masses of the light neutrinos. 

The paper is organized as follows. In Sec. Ull we introduce our formalism and set up the 
analytical framework under which the symmetry relations may be examined. In Sec. IIIII 
and Sec. HVt we test the texture zeroes of M^, and Mm respectively, numerically as well as 
analytically. In |V] and IVIl we examine the fi — t exchange symmetry and 5*3 permutation 
symmetry for M^, and Mm respectively. Sec. IVIII concludes. 
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II. THE ANALYTICAL FRAMEWORK 



A. parameterization of neutrino masses and mixing 

We parameterize the absolute values of neutrino masses in terms of three parameters mo, 
e and p as [3] 

mi = mo(l - p)(l - e) , ms = mo(l - p)(l + e) , m3 = mo(l + p), (8) 

where mg sets the overall mass scale of neutrinos, while the dimensionless parameters p and 
e can be expressed in terms of the solar and atmospheric mass scales as 

_ Am^tm _ Am| 

^ " ' ^ " 4m2(l - p)2 • 

Clearly, p is positive (negative) for normal (inverted) mass ordering of neutrinos. The sum 
of neutrino masses may be expressed in terms of the above parameters as 

5^m, = 3mo (l-0 ^ 1 eV. (10) 

i 

The condition < mj < 0.5 eV then yields 

mo > 0.025 eV , 2.43 x 10"^ < |p| < 1 , 8 x 10^^ < e < 1 . (11) 

The value of |p| approaches unity as mo approaches its lowest allowed value. The value of e 
can be > 0.01 only for normal mass ordering and mo < 0.06 eV, whereas e <^ |p| everywhere 
except for mo ~ 0.025 eV. Taking the best-fit values of solar and atmospheric neutrino 
masses, in Fig. [1] we show the values of p and e as functions of mo for normal as well as 
inverted hierarchies. For the purpose of this paper, we divide the neutrino parameter space 
into three scenarios: 

(i) Normal mass ordering with hierarchical masses (NH), where mi -C m2 -C m^. The 
current data give p ~ 0.85 and e ^ 0.92 in the extreme limit, however these values decrease 
rather rapidly as mo increases, as can be seen from Fig. [T] This scenario can then be analyzed 
through a perturbative expansion in the set of the small parameters p = 1 — p and e = 1 — e 
in the extreme limit, however one has to be careful while treating quantities like (1 — p^), 
which stays higher than 0.3. 

(ii) Inverted mass ordering with hierarchical masses (IH), such that ms <^ mi < m2. In this 
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FIG. 1: The parameters p and e as functions of mo for (a) normal ordering and (b) inverted 
ordering of neutrino masses, for best-fit values of Am^tm and AuIq. 



case, e ^ 1 and |l + p| ^ 1, so we can use a perturbative expansion in the small parameters 
p = 1 + p and e. 

(ill) Quasidegenerate neutrinos (QD), where mi ~ m2 — ms, with either mass ordering. In 
this case, |p|, e ^ 1, so that we can use these two quantities as small parameters. 

In addition, at appropriate places we shall also consider ^13 and 623 = 623 — it/ 4 as small 
parameters in order to facilitate a perturbative expansion. 



B. Elements in M^, M'^, and Mm 



Let the low energy neutrino mass matrix in the flavor basis be written as 



^ a h c"^ 

M^= b d e . (12) 
[c e f) 

Since the neutrinos are Majorana, Mj, is symmetric. In terms of the parameterization of 
neutrino masses in Sec. Ill At mixing angles and CP violating phases, the elements of Mi, 
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may be written as 



a — moe 



e2^^i(-l + p)(-l + e)4c 



13 



e-"'(l + p)s23Si3 

+e2^<^2(_i + + e)5i2 (-C12C23 + 6^-^512523513) 



e-''^(l + p)c23Si3 

_e2^^l(-H-p)(-l + e)Ci2 (-S12523 + e^'ci2C23Sl3) 

+e2^-^^(-l + p)(l + e)si2 (C12S23 + e''c23Si2Si3) 



(i + p)4s^ 



23 



+e2^'^H-l + P)(-l + e) (C23S12 + e''ci2S23Si3)' 



■1 + p)(l + e) (C12C23 - S12S23S13)' 



^^gi(XM+X.) 



(l + p)c23Ci3S23 

_(,2i<t>i(^_l ^ ^ g) ^g^^g^g _ e^''ci2C23Si3) (C23S12 + e^''ci2S23Sl3) 

^g2i02^_l + p)(l + e) (C12S23 + e*^C23Si2Si3) (C12C23 - e*'^Si2S23Si3) 



2«Xt 



(1 + p)c^3C?3 

+e2^'^i(-l + p)(-l + e) (S12S23 - e^'ci2C23Si3) 
-e^i^^{-l +p){l + e) (C12S23 + e''c23Si2Si3)' 

The inverse of the neutrino mass matrix, M~^, can be written as 



^ A B 



M 



Det{M^) Det{M^) 



B D E 
C E F 



where 



Det{M,) = mg(l - p)'(l + p)(l - e')e 



is the determinant of My, with = Xe + X^l + Xr- Here M is the adjoint neutrino mass 
matrix. From Eq. f |T4l) it is obvious that texture zeroes in M~^ are the same as those in M. 
The elements of M can be written in terms of the masses and the elements oiU = t/pMNS 
matrix as 



A 



B 



mae-^-'— ^ [(1 - pfil - e') (f/2if/32 - U22U^i] 



,2^2i(x„+Xr) 
„2^ 



c 



^2gi(xe+XM+2XT) 
+ (l-p')(l-e 

^2gi(xe+2x^+XT) 



+ (1-P')(1- 
+ (1-P')(1 + 



2„2j(x.+Xr) 



+ (l-p')(l+e 

^2gi(2xe+Xn+XT) 

+ (l-p')(l + e 



(f/2lf/33 - f/23f/3l)' + (1 - P')(l + e) (f/22f/33 " f/23t/32)'] 
(1 - pf{l - e^) (f/2lf/32 - f/22f/3l) (t/3lt/l2 " t/llt/32) 
(t/3lf/l3 - f/33f/ll) (f/2lt/33 " f/23f/3l) 
(t/22f/33 - f/23f/32) (f/32f/l3 - U^^)] 

(1 - p)2(l - e^) (f/2lf/32 - f/22f/3l) (t/22f/ll " U12U21) 
{U12U23 - f/l3f/22) (f/22f/33 " U^2U23)\ 



[(1 - pf{l - e^) (f/nf/32 - Ui2Uuf + (1 - p')(l - e) ([/iif/33 - t/i3f/3l)' 



([/i2f/33 - f/32?7i3)'] 

(1 - - e^) (f/nf/32 - f/l2f/3l) (f/l2f/21 - f/llt/22; 

(f/llf/33 - f/l3f/3l) (?72lf/l3 - f/llt/23) 
(f/l2f/33 - f/l3f/32) (f/22f/l3 - f/l2f/23)] 



2„2j(x.+XM) 



[(1 - p)2(l - e^) (f/nf/22 - t/i2f/2i)' + (1 - p')(l - e) (t/iit/23 - UMif 



+ (l-p2)(l+e) (f/i2f/23 - f/l3f/22)'] 



(16) 



In order to analyze the heavy majorana neutrino mas matrix Mm-, we invert Eq. ([7]) to 
obtain 

Mm = -mlM-^rriD • (17) 

Following 1^, we choose the "flavor" basis for heavy Majorana neutrinos in which the Dirac 
mass matrix is real and diagonal, 



niD = diag(x,|/, z) 



(18) 
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In this basis, Mm may be written as 



x^A xyB xzC 
xyB y^D yzE 



\ 



Mm 



1 



(19) 



Det{M^) 




I 



Again, the texture zeroes of Mm are same as the texture zeroes of M~^, and hence those 
of M, which have relatively tractable analytical expressions. The discrete symmetries like 
— r exchange or 5*3, on the other hand, also depend on the values of the Dirac masses. 
However even in that case, we can test for certain relations between elements of Mm that 
are independent of these Dirac masses, as we shall see in Sec. I VII 

C. Quantifying deviation from exact symmetry in the bottom-up approach 

In the traditional top-down approach, discrete flavor symmetries like /i — t exchange, 
Ss-permutation, etc. are assumed in the neutrino mass matrix, which predict the mixing 
parameters measured in the low energy neutrino oscillation data. The main purpose of these 
symmetries is to understand why the (1-3) family mixing is small, while the (2-3) family 
mixing is almost maximal and the (1-2) family mixing is large but not maximal. Traditionally 
these symmetries are employed to set f/13 = 0. The dynamical breaking of these symmetries 
then may predict a non-zero f/13, which is then compared with the experimental value. Since 
the current low energy neutrino oscillation data have large uncertainties, the data allow an 
enormous freedom to propose such discrete flavor symmetries in the top-down approach. 

However it is also crucial to examine, by starting with the available low energy neutrino 
oscillation data, the parameter space of neutrino mass matrix where such discrete flavor 
symmetries can be realized. We call it the bottom-up approach. In this approach, since the 
low energy data have intrinsic uncertainties, the symmetry relations can only be said to be 
satisfied approximately. One therefore needs to quantify when one may declare the relevant 
symmetry to be allowed. 

In M^, the magnitudes of all the elements are expected to be ~ tuq, as can be seen from 
Eqs. ( fT3l) . taking into account that the sine and cosine of 6^12, 6*23 ~ and also assuming 
that (pi ~ 0{1) in the absence of any symmetry principle. If an element \Miy{i,j) \ is <^ mo, it 
is either an accidental cancellation or the signature of a discrete symmetry at work. We take 
the position that for a sufficiently small value of ^, the observation \M^{i,j)\/mo < ^ would 
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indicate that the symmetry that would make M^{i,j) = is present. We choose ^ = 10~^, 
which is motivated by the accuracy to which the mixing angles are currently known. It 
also indicates the extent to which we tolerate the breaking of the discrete symmetries under 
consideration. In other words, when \M,^{i, j)\/mQ < 10^^, we consider M^{i,j) to be 
effectively zero. Thus, we declare a texture zero viable if 



where the minimization is over all the allowed values {3a) of the mixing parameters at a 
particular value of mo. If this condition is not satisfied, then the symmetry that would lead 
to Mi,{i,j) = is ruled out. 

Similarly, from Eq. (fT6|) one can see that the elements of M are expected to be ~ tuq in 
the absence of any cancellations. Hence if 



then we conclude that the symmetry that requires M{i,j) = is still allowed. In the case 
of the discrete symmetries like fi — r exchange or 5*3, certain ratios are expected to be equal 
to unity. Here, we demand that the deviation of such ratios from unity to be less than 10~^ 
for the symmetry to be acceptable. Note that the right hand side of Eqs. (!20|) and (12T!) can 
be changed to any small number of one's choice, depending on how much deviation from the 
exact symmetry one is willing to allow. Our numerical results cover the complete relevant 
range, so the required numbers can be read off from our figures. 

Note that our criteria give the necessary conditions for a particular symmetry to hold. 
Further considerations may disallow some of the symmetry relations that are permitted by 
conditions in Eqs. fl20l) and fl2T|) . 

III. TEXTURE ZEROES IN 
A. Individual zeroes in Mi, 

Whether a particular element in the neutrino mass matrix can potentially vanish 
can be checked analytically from Eqs. (fT3ll . To simplify the expressions, we define three 
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FIG. 2: The minima of certain M^{i, j)/mQ are shown against mo for (a) normal and (b) inverted 
ordering of neutrino masses. The minima of the remaining elements of My/rriQ are less than 0.001 
for both normal as well as inverted ordering. The neutrino mixing parameters are varied in their 
3(T allowed range in this and the subsequent figures. 



quantities 



Ci = 


(1 - e)e^'^' , 


(22) 






(23) 


Cs = 


Size'' . 


(24) 



Note that since e > 0, we have \C,i\ < IC2I in case of NH, while |Ci| ~ IC2I ~ 1 in the IH and 
QD scenarios where e ^ 1. 

The following observations may be made from the analytic expressions for a, 6, c, rf, e, /: 
(i) a = M,(l, 1) = mo(l - p)e'''^^ [^4, + C2sj, + (1 + p)Cf] ■ 

In NH, while the (1 — p) provides one suppression factor, the other suppression is provided 
by the terms involving Q, all of which have phases that can be adjusted so as to cause a 
cancellation among the three terms inside the square bracket. Thus, M^(l,l) can vanish 
for NH as can be seen from Fig. [2l Note that if ^^13 were extremely small so as to make the 
(^3 term negligible, the cancellation could not have been achieved since current data implies 
IC1C12I < |C2"Si2| in this extreme hierarchical limit. So a significantly non-vanishing ^13 is 
crucial for allowing a texture zero of M,y{l, 1) for NH. 

In IH and QD, since |Ci| ~ IC2I ~ 1 and ~ -512 > 0.26, the minimum magnitude 
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of the sum of the first two terms in the square bracket is ~ 0.26. On the other hand, 
iCll < 0.04, so this term cannot cancel the first two. Thus, the smallness of 6*13, combined 
with the non-maximality of 9\2i ensures that My{l, 1) cannot be a texture zero for IH and 
QD scenarios. 

In the normal mass ordering, extreme hierarchy allows the texture zero of M,^(l, 1) 
whereas quasidegeneracy prevents it. The transition between these two extremes as a func- 
tion of mo may be observed in Fig. [2l 

(ii) h = M,(l, 2) ^ -(mo/v^)e^(^^+>^'')[^]f, + ^^23)] , 

c = M,(l, 3) ^ +(mo/y2)e*(^'=+^^) [n, + 0(0?^, ^^23)] , 

where 

n, = n, = {l- p)[Ci2Si2(Cl - C2) - C3(ClC?2 + C2^5?2)] + (1 + P)Q ■ 

In NH, if = then VL^ and Vt^ cannot vanish, since |Ci| < |C2|- However, with a 
non-zero (^3, the phases 0i,02 and 5 may be chosen properly to make VLb and Vic vanish. 
Note that the (1 + p) coefficient of the Q term makes that term comparable to the terms 
involving Ci,C2, which are suppressed by the coefficient (1 — p) for NH. Thus, My(l,2) and 
Mj,(l, 3) can be texture zeroes in NH as long as ^13 is not extremely small. 

In IH and QD, \C,i\ ~ IC2I, so that 0i ~ 02 can make Ci ~ C2 ~ 0. Thus even with ^13 
vanishing, M^(l, 2) and M^{1, 3) can be texture zeroes of the neutrino mass matrix. 

(iii) d = M,(2, 2) ^ {mo/2)e^^^^ [Q^ + 0{el„ 6,,)] , 

/ = M,(3,3) ^ {mo/2)e^'^^[Qf + 0(^?3>^^23)] , 

where 

nd = nj = ii + p) + ii- p)[(Ci4 + C2CI2) + 2C3Si2Ci2(Ci - C2)] . 

In NH, (1 + p) is as high as 1.9. This cannot be cancelled by the other terms involving 
Ci, C2 since these terms are already suppressed by (1 — p). As a result, M,^(2, 2) and M^{3, 3) 
cannot be texture zeroes for NH. This can be seen from Fig. [2|^a). In IH, (1 + p) ^ 0. 
While the coefficient of (1 — p) in and Qf cannot vanish for ^13 = 0, it can be made 
to vanish for ^13 7^ with proper choices of (f)i,4>2 and 6. On the other hand in QD, the 
choice of 0i = 02 = 7r/2 would make Qd and Qf vanish even when 6'i3 = 0. Thus Mj^(2,2) 
and Mj,(3, 3) can be texture zeroes in IH and QD scenarios, the former scenario requiring a 
nonzero ^13. 

(iv) e ^ M.(2, 3) ^ (mo/2)e^fe'+>^^)[(l + p) - (1 - p)(CiS?2 + (20^2) + 0{9l„ O^,)] . 
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TABLE I: Viability of individual texture zeroes for M^. Here indicates that the texture is 
allowed even for ^13 vanishing, Q indicates that the texture is allowed but needs a nonzero ^13, 
whereas (8) indicates that the texture is not allowed. 

In NH, arguments similar to above can be used to show that the large value of (1 + p) 
cannot be cancelled by the other terms that are suppressed by (1 — p). As a result, M^{2, 3) 
cannot be a texture zero for NH. In IH, while (1 — p) vanishes, the coefficient of (1 — p) term 
cannot vanish owing to the non-maximal nature of 612, preventing M^{2,3) from being a 
texture zero for IH. In QD, however, 0i = 02 = 7r/2 would make e vanish, so that M^{2, 3) 
can be a texture zero. 

In Table [T] we summarize our results on the viability of the texture zeroes in various 
scenarios. 



B. Two zeroes in M, 



In order to quantify the two-zero textures in in the bottom-up approach, we define 

■v/|M,(z,j)P + |M.(A;,/)|2 



M{ij -kl)= Min 



(25) 



There are fifteen M{ij — kl) constructed out of six Mi,{i,j) elements. From the discussion 
in Sec. IIII Al it is clear that the five M{ij — kl) involving M,y{l, 1) cannot vanish unless 
the mass ordering is normal and mo ^ 0.035 eV. Moreover, even in such a case the other 
element of the pair cannot be M,^(2, 2), M^(2, 3) or M;^(3,3). The only two-zero textures 
including vanishing M^{1, 1) are therefore 



X 
X X 
XXX 



and 



0X0 
XXX 
X X 



(26) 



where X denotes an element that may or may not vanish. That both these patterns corre- 
spond to NH can be seen from Fig. [31 The above two- zero textures with normal ordering 
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FIG. 3: M(ll — 12)/mo and M(ll — 13)/mo as functions niQ in (a) normal and (b) inverted ordering 
of neutrino masses. 



imply that the neutrinoless double beta decay (Oz//3/5), whose rate is proportional to 



(27) 



|M.(1,1)P= 5^m,[/i2 

r ' 

is unobservable since |M^(1,1)| < 10 ^mo [17|]. In other words, if the neutrino masses are 
normally ordered and the 0z//3/3 is observed then the corresponding symmetry which drives 
the cancellation in the neutrino mass matrix leading to these textures should be forbidden. 
The five textures 



X 
X X 
X X 



X X 
0X0 
X X 



XXX 
X 
X X 



XXX 
X X 
X 



(28) 



X X 
X X 
XJ 

cannot satisfy the neutrino data since the first two lead to 612 = 0, the second and the 
third lead to 6*23 = 0, while the last two are inconsistent with a small 613 and large 623 
simultaneously. 

The remaining five pairs of M^{i,j) textures are allowed by both normal and inverted 
ordering of neutrino mass pattern at large mo, i.e. they are allowed in the QD scenario. 
These textures are 



X X 
X 
XXX 



X X 
X X 
X X 



X X 
X X 
X X 



X X 
XXX 
0X0 



X X 
X X 



(29) 
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FIG. 4: M(ij — kl)/mQ as functions of tuq for (a) normal and (b) inverted ordering of neutrino 
masses. The fluctuations seen in this and some of the subsequent figures are numerical artifacts, 
a consequence of the inability of the randomly chosen neutrino parameters to find the actual 
minimum. 

The mo values at which these textures start becoming viable can be read off from Fig. |H 
Notice that the seven allowed two- zero textures, given by Eqs. ( !26|) and (!29l) . have already 



been predicted in the top-down scenario 



Here in addition to obtaining them through a 



bottom-up approach, we have also correlated them with the measured values of tjiq and 6*13. 



C. Three zeroes in M,, 



In order to quantify the possible three-zero textures in M^, in the bottom-up approach, 
we define 

V|M,(z,j)P + |M,(fc,/)|2 + |M,(m,n)|2 



M(ij — kl ~ mn) = Min 



(30) 
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From the two-zero textures obtained in the last subsection, one can deduce that the only 
possible three-zero textures of M^, are the combinations (12 — 22 — 33) and (13 — 22 — 33). In 
Fig. O we show these two combinations M{ij — kl — mn)/mQ as functions of mo for normal 
as well as inverted ordering of neutrino masses. It can be seen that M{ij — kl — mn)/mQ 
does not vanish in any region of the neutrino parameter space. Thus, there are no viable 
three-zero textures for M,y. 
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FIG. 5: M{ij — kl — mn)/mQ as functions of rriQ for selected values of {ij,kl,mn) for (a) normal 
and (b) inverted ordering of neutrino masses. For all the rest (ij, kl, mn) values, M{ij — kl — mn) 
is always greater than O.lmo and hence the corresponding symmetry is forbidden. 



IV. TEXTURE ZEROES IN M~\M AND Mm 

The texture zeroes in M~^, M and Mm are identical, as can be seen from Eqs. f|T4|) and 
(|T9l) . The quantification of "smallness" of the matrix elements is even more arbitrary for Mm 
due to the presence of the unknown Dirac masses x, y, z. However for the sake of uniformity, 
we determine the texture zeroes of M through the quantitative criteria described in Sec. Ill Ci 
and apply the same criteria for the zeroes of Mm- This correspondence preserves all the 
zeroes in M, and does not add any additional ones due to the hierarchy of Dirac masses, for 
example. 

We thus continue with determining the zeroes of M on lines similar to the last section 
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where we determined zeroes of M^. For convenience, we define 

M,., - « - . M,„ ,\/|g('.^)l-^|g(M)|-.|ffK.)|A ^ 

The quantities Min(|M(?, M (ij — fc/) and M{ij — kl — mn) need to be less than 10~^mo 
in order to quahfy as one-, two- or three- zero textures respectively. 



A. Individual zeroes in M 
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FIG. 6: Individual elements of Min(M(i, j))/mQ as functions of mo for (a) normal and (b) inverted 
ordering of neutrino mass spectrum. The minima of the remaining elements which are not shown 
in the figure are less than 0.001 for all values of mo- 



(i) A ^ M(l, 1) ^ my^^^^+^^) [(1 - p'){Cisl, + C2C?2) + (1 - p)'CiC2C|] • 

In NH, the last term is suppressed quadratically in (1 — p) as well as in (^3, whereas the 
(1 — p^) term does not have such a strong suppression. Moreover, the coefficient of (1 — p^) 
cannot vanish since IC1C12I IC2Si2l- As a result, M (1, 1) does not vanish for NH. In IH, 
p^ ~ 1, so that the leading term vanishes. Further, the (1 — p)^ dependent term vanishes at 
^13 = 0. M(l, 1) as a texture zero is therefore allowed. In QD on the other hand, the (1 — p^) 
term and its coefficient are nonzero, while the second term is suppressed by (1 — pYOf^ and 
is unable to cancel the first term. As a result, M(l, 1) cannot vanish for QD. 
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In the inverted mass ordering, the hierarchical hmit allows a texture zero of M(l, 1) 
whereas the quasidegenerate limit prevents it. The transition between these two limits as a 
function of tuq can be seen in Fig. [6](b) . 

(ii) B = M(l, 2) ^ +K/v/2)e*fe+>^''+2x.)[^^ + ^(^2^^ ^ 
C = M(l,3) ^ -(m2/v/2)e*(>^'=+2xM+x.)[^]^ + OiOl^O^s)] , 

where 

= ^^C = (1 - P')[C12S12(C1 - C2) + C3(ClsL + C24)] - (1 - p)'ClC2C3 . 

In NH, the last term is suppressed by (1 — p)'^6is, and is therefore rather insignificant. 
However the coefficient of the (1 — p^) term cannot vanish with any choice of the Majorana 
and Dirac phases, so M,y{l,2) and M^(l,3) cannot be texture zeroes in NH. This feature 
can be seen from Fig. [6](a). In IH and QD, since |Ci| ~ IC2I5 one may choose 0i ~ 02 to 
make Ci ~ C2 ~ 0, so that M(l, 2) and M(l, 3) can be texture zeroes. 

(iii) D ^ M(2, 2) ^ my^(^^^^^)[nn + 0(61,, §1,)] , 

F = M(2, 3) ^ m2e2^fc+XM)[^]^ + 0{ef^, ^fg)] , 

E = M(2, 3) ^ m2e^(2x.+XM+xr)[^g + (9(^2^^ ^2^)] ^ 

where 

l^D = ^^^ClC2 + ^^^(Cl4 + C25L - 2C3(Cl - C2)C12S12) , 

= ^^^ClC2 + ^^^(ClcL + C2SI2 + 2C3(Cl - C2)C12S12) , 

-2(1 - p)'ci2Si2ClC2C3 
= -^^^ClC2 + ^^^(ClC?2 + C2SI2 + 2C3(Cl - C2)C12S12) • 

In NH, the first term in each of the expressions Qd, ^f, is suppressed by (1 — p)^Ci) 
whereas the coefficient of (1 — p^) can be made to vanish with an appropriate choice of 
phases, as long as 61^ is finite and Cs can participate in the cancellation. So M(2, 2), M(2, 3) 
and M(3, 3) can be texture zeroes for NH. In IH, though the (1 — p^) term vanishes, the 
coefficient of the (1 — p)^ term is of the order of unity and as a result, none of the above 
three can be texture zeroes. In QD, the first term is suppressed by (1 — pYCii while the 
coefficient of (1 — p^) becomes 

(ClC?2 + C2S?2 ± C1C2) 
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A B C D F E 
NH ® ® ® Q Q Q 

IH yj ® ® ® 

QD ® V V V V V 
TABLE II: Viability of individual texture zeroes for Mm- Here ^J indicates that the texture is 
allowed even for ^13 vanishing, Q indicates that the texture is allowed but needs a nonzero ^13, 
whereas ® indicates that the texture is not allowed. 

when ^13 = 0. The choices of the Majorana phases 0i = 02 = 7r/2 and 0i = 02 = make 
the above expression vanish for the + and - sign respectively. Thus, Mj,(2, 2), Mj,(2, 3) and 
M^(3, 3) are allowed as texture zeroes for QD. 

In the inverted mass ordering, the hierarchical limit allows a texture zero of the above 
three matrix elements, whereas the quasidegenerate limit prevents it. The transition between 
these two limits can be seen in Fig. EJ^b). 

In Table [IT] we summarize our results on the viability of the texture zeroes in various 
scenarios. 



B. Two zeroes in M 



Out of fifteen elements in M{ij — kl), five two-zero textures involve a zero of M(l, 1). 
These are clearly not allowed for NH and QD, since the value of M(l, 1) itself is high in these 
scenarios. In IH, M(l,l) can be small, however M (2, 2), M(2, 3) and M(3,3) are large, so 
that M(ll — 22), M(ll — 23) and M (11 — 33) cannot be texture zeroes. This leaves us with 
the two possible textures 



X 
X X 
XXX 



0X0 
XXX 
X X 



(33) 



However, these two inverse neutrino mass matrices correspond to the last two two-zero 
textures of M^, matrices in fl28l) . which are forbidden. Thus, even the textures in fl33l) are 
not allowed. 
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FIG. 7: Minima of M{22 — 23) and M(23 — 33) as functions of mo for (a) normal and (b) inverted 
ordering of neutrino masses. 



X 
X X 
X X 



^ ^ X X 0^ 



J 



X X 
OCX 



(34) 



The textures 

X X 
0X0 
yX 

are ruled out since in each of these, one of the neutrinos is completely decoupled, whereas 
we need large values for 6^12 as well as 6*23. 

Textures involving zeroes of M(2, 2), M(2, 3) and M(3,3) can be allowed only for NH 
and QD. It is observed that M(22 — 23) and M(23 — 33) can vanish in the NH scenario, but 
not in the rest. The dependence of the vanishing of these two textures. 



^ X X x\ fx X x^ 

X 
X X 



X X 
X 



(35) 



can be seen from Fig. [7| where we have plotted M{ij — kl)/mQ as a function of mg for normal 
as well as inverted ordering of neutrino masses. 
The remaining five two-zero textures. 



X X 
X 
XXX 



X X 
X X 
X X 



^ '^X X 0^ 



X X 
X X 



'^X X 0^ 



XXX 
0X0 



and 



XXX 
X X 
X X 



\ 



(36) 
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FIG. 8: Selected M{ij — kl)/m'^ as functions of mo for (a) normal and (b) inverted ordering of 
neutrino masses. 



are allowed by both normal and inverted ordering of neutrino mass pattern, though only in 
the quasi-degenerate limit. This can be seen quantitatively in Fig. [HI Note that the textures 
given in Eqs. fl5Bl) . and fl5B]) are already predicted in the top-down scenario Here in 
addition to obtaining them through a bottom-up approach, we have also correlated them 
with the measured values of mo and 6*13. 



C. Three zeroes in M 



Using the two-zero textures of M predicted in the last section, there are only two possible 
three-zero textures that are viable: (12 — 22 — 33) and (13 — 22 — 33). Fig. [9]shows the value 
of M{ij — kl — mn)/m^ for these two cases. As can be seen from the figure, there are no 
three-zero textures of the adjoint neutrino mass matrix M, and hence of Mm- 



21 



(a) (b) 
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FIG. 9: M{ij — kl — mn)/mQ as functions of mo, for potentially viable values of [ij, kl, mn) in case 
of (a) normal as well as (b) inverted ordering of neutrino masses. 



V. /i - r AND 5'3 SYMMETRIES IN My 
A. fJ- ^ T exchange symmetry in Mi, 

A fi — T symmetry in M^, implies that, in the notation of (IT^ . we have d = f and b = c. 
The deviations from fi — t symmetry may be then quantified in terms of the dimensionless 
parameter 



A 



flT 



d 


2 


1-^ 




+ 








C 



(37) 



The existence of a — r symmetry in Mj,, as per our convention in Sec. Ill C\ requires 
Min(A^^) < 10~^. In order to analytically understand the constraints on various parameters, 
we expand A^^ in terms of the relevant set of small parameters. 
For NH, using Eq. (|T3|) one can calculate 



;i - (i//)NH 

(1 - 6/c)nh 



I _ g«(XM-XT) 



2^13 + e^(2<^2+'^)psin(2^i2^ 



2^13 -e^(2</'2+5)psin(20i2^ 



23 J 



(38) 



For A 



0, we need both the above expressions to vanish, which can only happen if 
IT and 6*13 ^ psin(26'i2). 
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In the IH scenario, we get 

(l-6/c)iH = l + e^(>^--^^) + 0(p,e,^i3,^: 



23 J 



(39) 



so that \Xft ~ Xt\ ~ tt is enough to get A^^ pa to zeroth order in the relevant small 
parameter. 
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FIG. 10: The contours of A^,-, the deviation from the /U — r exchange symmetry in M^^ in the 
plane of sin 6*23 versus sin for (a) normal and (b) inverted mass ordering, for mo = 0.03 eV. 



The QD scenario also gives 

(1 - c?//)qd 
(1 - Vc)qd 



(40) 



so that we expect A^^ ^ to be satisfied if — Xt| ~ tt. 

In Fig. [TOl we show the numerical results for Min(A^T-) for hierarchical neutrino masses 
with both mass orderings. The current data allow for exact fi — t symmetry whatever the 
value of mo or nature of mass ordering, it will be possible to rule out this symmetry if a 
significant value of ^13 is measured at experiments. 



B. S3 permutation symmetry in Mj 



The 5*3 permutation symmetry in demands that, in addition to the conditions d = f 
and b = c satisfied by the [i — t exchange symmetry, one also needs h = e and a = f. The 
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deviation from 5*3 symmetry then can be quantified by the dimensionless parameter 



d 
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b 
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a 
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1-^ 
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+ 




+ 






C 




'-J 




e 



(41) 



Clearly, the jj, — t symmetry is a subset of the 6*3 permutation symmetry. Then the 
discussion in Sec. IV Al shows that we need the conditions in Eqs. (1551) . fl5^ and fHUl) satisfied 
in the NH, IH and QD scenarios respectively, to ensure |1 — d/ f \ <^ 1 and |1 — h/c\ ^ 1. 
In addition, we need to check if the terms |1 — a//| and |1 —h/e\ can be small quantities. 
We shall again analyze the relevant parameter space by considering the analytic expansion 
in the relevant small parameters as done in the previous section. 

In the NH scenario, 

(l-a//)NH = l + 0(p, 6,^13, ^23) , 

(1 - 6/e)NH = 1 + 0{p\ e^ el^, el,) . (42) 

None of the above conditions can be satisfied in this region of parameter space, where p -C 1. 
Therefore, one cannot have 5*3 symmetry. 
In the IH scenario, 

/ 2j0i 2 I „2j02o2 \ 

(1 - a/fU = 1 - 2e-(---) + Oip, e, ^,3, ^^23) • (43) 

(/ 2j0i _ 2102^ e \ 
,2.,.,2^ + e4"^" J + Oip, e, ^,3, ^~23) . (44) 

These conditions correspond to strict correlations between the flavor phases, the Majorana 
phases and the mixing angle 612- However, the above two equations cannot be satisfied 
simultaneously. This can be seen as follows: the simultaneous conditions a// ~ 1 and 
6/e ~ 1 would imply {p/eY ~ a//, or {hjef — a/ f ^ However, one gets 

[{h/ef - a/fU = (^..,,^2^^,2.,,,2^)2 + ^(P' ^13, ^23) , (45) 

which does not vanish even to the zeroth order of the small parameters. Thus, the current 
data do not allow S3 symmetry to hold in the IH scenario. 

On the other hand, when neutrinos are quasi-degenerate in mass, 

(1 - a/f),n = 1 - 2e-(---) {j^^^^^^^J + Oip, e, 9,,, ^,3) ■ (46) 

1 _ ,2.,.,2^ _ ^2.^2 J + 0{P, 6, ^,3, ^^23) • (47) 
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FIG. 11: The deviation of S3 symmetry in M^^ as a function of mg for (a) normal and (b) inverted 
mass ordering. 
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which also constrain 612, 0i, 02, Xe, Xr- The above two equations have simultaneous solutions, 
so the current data allow 5*3 symmetry to hold in the QD scenario. 

Thus, the mo value determines whether 6*3 is an allowed symmetry in the normal as well as 
inverted ordering. Fig. [TT] shows the transition between hierarchical and degenerate scenarios 
where 5*3 becomes viable. In Fig. [T21 we show the deviation from 6*3 symmetry, which we 
have quantified in terms of A53, as a function of 612 and (pi. Since these two parameters can 
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be directly constrained from the neutrinoless double beta decay, these experiments will be 
crucial for testing the 5*3 symmetry in neutrinos. 



VI. fi-T AND S3 SYMMETRIES IN Mm 
A. I-L — T symmetry in Mm 

A — r exchange symmetry in Mm implies, from fll9p . that xyB = xzC and y^D = z^F. 
Since we are completely ignorant about the Dirac masses x, z, we cannot use these two 
conditions separately. However, we can use the Dirac-mass independent combination of 
these conditions, which gives {B/CY = (D/F). The deviation from this fi — t symmetry 
relation may be quantified by 



Since this condition is less restrictive than that for the ^-t symmetry in M^, which itself 
is consistent with the current data irrespective of the value of mo, one would expect that 
the /i-r symmetry will hold also for M. Indeed, it is found that Min(A^^) ^ 10^^ for all 
allowed values of the neutrino parameters. This can be understood as follows. In the NH 
scenario, the expansion in terms of the small parameters ^^13, ^^23, e, P gives 

{B'/C^ - D/F)^n = 0{~p\ e^ el,, ~el,) , (49) 

where we have used Eq. (fT6|) . Note that even terms linear in the small parameters 6'i3, 6^23, e, p 
are absent. This indicates that in this scenario, the p — r exchange symmetry is easily 
satisfied. In the IH scenario, 

which vanishes for 6'23, Q\z — 0. In the QD scenario also, p — t exchange symmetry is easily 
satisfied since 

(B^C^ - D/F)qu = 0{p, e, ^13, ^~23) . (51) 

Thus the p — r symmetry is consistent with the current data irrespective of rriQ or the mass 
ordering of neutrinos. 

In the IH scenario, the deviation from the p-r symmetry would manifest itself in the 
angles 613 and 6'23. Clearly the symmetry is valid if ^13 = = 6*23. If ^^13 = exactly and 
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FIG. 13: The contours of A^,-, the deviation from the fi-T permutation symmetry in Mm, in the 
plane of sin 6*23 versus sin 6*13 for (a) normal and (b) inverted mass ordering, for mo = 0.025 eV. 



6*23 7^ 0, Eq. (I5U|1 indicates that — D/F)iii cannot vanish, thus the symmetry is 

not obeyed. However, if the upper bound on 613 is significantly nonzero, the two terms in 
Eq. fl50l) can cancel each other since the ^13 term is enhanced by the factor (1/p). Thus to 
satisfy /i — r symmetry in Mm, either ^^13 ^ ^ ^^23 7 or the upper bound on 6*13 should be 
large. In the latter case, the actual value of sin^^i3 would have to be much smaller than p to 
prevent the second term from becoming too large. These features can be seen from Fig. [T31 
Indeed, for sin6'i3 0.01, it is possible for Mm to satisfy the fi-r symmetry for any value of 
6^23. For smaller values of 6^13, one needs 6*23 ^ 1.5°. 



B. 53 symmetry in Mm 

The 5*3 symmetry in Mm implies, from f|T9l) . that xyB = xzC = yzE, x^A = y'^D = z^F. 
The Dirac-independent conditions that can be obtained are (B/C)^ = (D/F) as in the case 
of the fi — T exchange symmetry, and the additional condition [B / E^ = {A/ F) is required 
for 5*3 symmetry to be satisfied in Mm- Therefore, the deviation from 6*3 symmetry can be 
quantified as: 

A53 = 

In the NH scenario, 
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(52) 



so that, combined with Eq. (H9i) . it is conceivable that 5*3 symmetry holds. 
In the IH scenario, we get 



{B''/E''-A/FU = 0{p,e,d,M ■ 



(54) 



The extent to which 5*3 is satisfied then depends only upon the extent to which the /i-r 
symmetry is obeyed. The effect of the deviation from this symmetry on the mixing angles 
is exactly the same as given in the /i-r case. The 5*3 permutation symmetry in Mm can 
thus be satisfied for either sin6'i3 >, 0.01, or |6'23| ^ 1.5° and sin^^i3 <, 0.01, as can be seen 
in Fig. [131 In the latter case, the actual value of sin ^13 needs to be much smaller than p. 
For degenerate neutrino masses, 



Since at the moment we have no restrictions on the Majorana phases, proper choice of 
the values of these phases would allow us to make {B'^/E'^ — A/F)qd vanish. Combining 
this with Eq. fIFIl) . one can claim that it is possible to have A53 ^ for quasi- degenerate 
neutrinos. 

VII. CONCLUSIONS AND OUTLOOK 

We examine texture zeroes and discrete fiavor symmetries in the neutrino mass matrix 
through a bottom-up approach. We develop a formalism that uses the low energy data, 
checks its consistency with the desired texture or symmetry, and quantifies the deviation 
from this symmetry. To this end, we parameterize the neutrino mass matrix in terms of 
the masses, mixing angles, CP violating Dirac and Majorana phases, as well as the "fiavor" 
phases that have no phenomenological implication but may be fixed by the mechanism of 
neutrino mass generation. We consider both the normal as well as inverted mass ordering 
of neutrinos, and the cases when the neutrino masses are hierarchical and quasidegenerate. 
The results can be described analytically in the three scenarios: normal ordering with hi- 
erarchical neutrino masses (NH), inverted ordering with hierarchical neutrino masses (IH), 
and any mass ordering, but quasidegenerate neutrinos (QD). In each independent scenario, 
we identify parameters that are small and hence can be used in a perturbative expansion. 




(55) 



g2j</)ig2 _^ e2i</-2s2 _ g2i(0i+02) 
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This provides a universal framework to analytically understand many results that were only 
known numerically before. In addition, it allows us to make new predictions about symme- 
tries that can be verified by numerical means. 

In the bottom-up approach, the extent to which the symmetries can be tested is limited 
by the accuracy of the experimental data on neutrino parameters, and hence one can only 
talk about the symmetry being approximately satisfied, or it being consistent with the data. 
In order to achieve this, we quantify the deviation from such a symmetry in each case 
through a quantity that vanishes in the limit of exact symmetry. The minimum value of 
such a quantity that is consistent with the current data is an indication of the extent to 
which the symmetry is valid. If this minimum value is zero, the symmetry is clearly obeyed. 
However, our formalism allows even for quantification of the breaking of the symmetry in a 
model independent, unified way. For illustration, in each scenario we consider the symmetry 
to be viable if the relevant quantity is less than 10~^. 

It is found that the viability of texture zeroes strongly depends on the absolute neutrino 
mass scale, the nature of mass ordering, and the angle ^13. For the neutrino mass matrix 
M^, all the six possible one- zero textures are allowed, though only in certain scenarios. For 
example, in normal mass ordering one can have only three possible one-zero textures when 
the neutrino masses are hierarchical, and five one-zero textures if they are quasidegenerate. 
The inverted ordering allows for four one-zero textures in the hierarchical limit (IH) and five 
in the quasidegenerate limit. We find that seven two- zero textures are allowed in M^, while 
no three-zero texture is permitted. 

The texture zeroes of the inverse mass matrix correspond to those of the Majorana 
mass matrix Mm of the heavy neutrinos in the Type-I seesaw mechanism, so we can treat 
them together. We find that all six one-zero textures are allowed for Mm, though only 
five of them are allowed for normal ordering. For inverted ordering, only three are allowed 
if the neutrinos have hierarchical masses, while five are allowed when the neutrinos are 
quasidegenerate. Seven two-zero textures are found to be consistent with the current data, 
while no three-zero texture is permitted. 

All the results with texture zeroes can be understood analytically in terms of the expan- 
sion in small parameters that we have introduced in each scenario. In most of the scenarios, 
the allowed freedom in choosing Dirac and Majorana phases helps us to take the mass ma- 
trix closer to the exact symmetry. Wherever one cannot do that, there is a special reason 
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that can be understood in terms of our analytic framework. Since the viabihty of the tex- 
ture zeroes is very sensitive to the absolute mass scale mo and the mixing angle ^13, future 
experiments aimed at determining these quantities will have a large impact, which can be 
gauged from our results. 

We use the same technique to examine flavor symmetries like /i-r exchange and S3. We 
define quantities involving elements of the M^, and Mm matrices that vanish identically for 
the exact symmetry. In order for the symmetry to be viable, we require this quantity to 
vanish to the zeroth order of the relevant small parameters. 

In Ml,, we find that the n — t symmetry is viable in IH and QD scenarios as long as the 
unphysical phases are related by \x^l ~ xA ~ ^r. In the NH scenario however, one needs 
\Xn ~ Xt\ ~ ti" and 613 <^ (1 — Aml^^/iml). As far as the symmetry is concerned, the 
current data allows it only if the neutrinos are quasidegenerate. 

Since the Dirac masses of neutrinos in Type-I seesaw mechanism are unknown, the con- 
straints put by the discrete fiavor symmetries on Mm are slightly weaker. As a result, the 
IJ> — T exchange symmetry is viable for NH and QD, while in the IH scenario one needs 
^13 -C |1 + Aml^j^J{Aml)\. The 5*3 symmetry in Mm also holds under these conditions. 

This paper provides a universal formalism to analyze the symmetries of neutrino mass 
matrices and illustrates its applications to a set of symmetries, viz. texture zeroes and 
discrete fiavor symmetries. Currently many of the texture zeroes and fiavor symmetries 
considered here seem to be viable. However, improved data on the mixing angles and masses 
expected in near future will be able to rule out some of these. Moreover, the formalism 
developed here can be apphed to a much wider class of problems related to the structure of 
neutrino mass matrix, that can give us more insights into the mechanism of neutrino mass 
generation. 
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